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ABSTRACT 

We present the largest-scale comparison to date between observed extragalactic X-ray binary (XRB) 
populations and theoretical models of their production. We construct observational X-ray luminosity 
functions (oXLFs) using Chandra observations of 12 late- type galaxies from the Spitzer Infrared 
Nearby Galaxy Survey (SINGS). For each galaxy, we obtain theoretical XLFs (tXLFs) by combining 
XRB synthetic models, constructed with the population synthesis code StarTrack, with observational 
star formation histories (SFHs). We identify highest-likelihood models both for individual galaxies 
and globally, averaged over the full galaxy sample. Individual tXLFs successfully reproduce about 
half of oXLFs, but for some galaxies we are unable to find underlying source populations, indicating 
that galaxy SFHs and metallicities are not well matched and/or XRB modeling requires calibration 
on larger observational samples. Given these limitations, we find that best models are consistent 
with a product of common envelope ejection efficiency and central donor concentration ~ 0.1, and 
a 50% uniform - 50% “twins” initial mass-ratio distribution. We present and discuss constituent 
subpopulations of tXLFs according to donor, accretor and stellar population characteristics. The 
galaxy-wide X-ray luminosity due to low-mass and high-mass XRBs, estimated via our best global 
model tXLF, follows the general trend expected from the Lx - star formation rate and Lx - stellar 
mass relations of Lehmer et al. (2010). Our best models are also in agreement with modeling of the 
evolution both of XRBs over cosmic time and of the galaxy X-ray luminosity with redshift. 
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1. INTRODUCTION 

Binary stars constitute a substantial fraction of stellar 
populations (SPs). In galactic fields and low density ag- 
glomerations such as open clusters between ~ 40% and 
~ 75% of stars are binaries (Duquennoy & Mayor 1991; 
Fischer & Marcy 1992; Fan et al. 1996; Raghavan et al. 
2010; Sana et al. 2012). In fact, stellar binarity may well 
be a universal characteristic of stellar evolution, since 
many single stars may either have been through a bi- 
nary phase, only to be ejected later, or be the result of 
a binary merger (Sana et al. 2012). Compared to sin- 
gle stars, binary stars are hosts to a range of additional 
processes (mass and angular momentum transfer, wind 
accretion, Roche-lobe overflow, common envelope ejec- 
tion etc.), making them ideal astrophysical laboratories 
for a whole range of physics not represented among sin- 
gle stars. Some of the most interesting processes that can 
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be probed are associated with accretion onto primaries 
that are compact objects (neutron stars and black holes). 
Due to the extreme energies involved, such cases are ob- 
servationally identified as X-ray binaries (XRBs). 

There have been many observational studies of XRB 
populations in external galaxies. Earlier work with 
the Einstein satellite showed that the X-ray emission 
is dominated by XRBs with low-mass (LMXBs) or 
high-mass donors (HMXBs) in early and late-type 
galaxies, respectively (Kim et al. 1992). Using Chan- 
dra's sub-arcsecond resolution, this result has now 
been established for individually detected XRBs in 
nearby galaxies (e.g. Kong et al. 2002; Soria & Kong 
2002; Trudolyubov et al. 2002; Sivakoff et al. 2003; 
Kim & Fabbiano 2003; Gilfanov et al. 2004; Gilfanov 
2004; Kim & Fabbiano 2004; Zhang et al. 2012; 
Binder et al. 2012, 2013, see also Fabbiano (2006) 
and references therein). Due to the longer evolutionary 
timescales of LMXBs, their integrated X-ray emission is 
an indicator of a galaxy’s total stellar mass (Gilfanov 
2004; Bogdan & Gilfanov 2010; Zhang et al. 2011; 
Boroson et al. 2011). In contrast, relatively short-lived 
HMXBs probe a galaxy’s star formation rate ((SFR), 
Grimm et al. 2003; Ranalli et al. 2003; Gilfanov et al. 
2004; Persic & Rephaeli 2007; Shtykovskiy & Gilfanov 
2007; Lehmer et al. 2010; Mineo et al. 2012). These 
observations have also established a break in the XLF 
of LMXBs at low luminosities {Lx ~ 10 37 erg s _1 , 
Gilfanov 2004; Revnivtsev et al. 2008, 2011) both in the 
Milky Way and in external galaxies. 

In the past, semi- analytical theoretical models have 
been introduced for the study of XRB populations. 
White & Ghosh (1998) and Ghosh & White (2001), as- 
sumed a time-dependent SFR and a simple rate model 
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to study the evolution of an arbitrary XRB population. 
These models predicted that the time required for bina- 
ries to reach the X-ray phase leads to a significant time 
delay between a star-formation episode and the produc- 
tion of X-ray emission from X-ray binaries from the pop- 
ulation. Wu (2001) created a simple birth-death model, 
in which the lifetimes of the binaries are inversely pro- 
portional to their X-ray luminosity, and calculated the 
XLFs of spiral galaxies. His models reproduce some 
features, such as the luminosity break in the observed 
XLFs of spiral galaxies. Piro & Bildsten (2002) argued 
that the majority of LMXBs in the field of elliptical 
galaxies have red giant donors feeding a thermally un- 
stable disk and stay in this transient phase for at least 
75% of their life. Most recently, in a series of papers, 
Bhadkamkar & Ghosh (2012, 2013a, b) started from stan- 
dard distributions of the parameters of those primordial 
binaries which are the progenitors of XRBs, and followed 
the transformation of these distributions with the aid of 
a Jacobian formalism as the binaries progress through 
different evolutionary phases. Following this methodol- 
ogy, they were able to derive estimates for the XLF and 
other population statistical properties, for both HMXBs 
and LMXBs. 

Binary population synthesis (PS) modeling codes can 
provide a unique tool to understand the physical prop- 
erties of this important population in a statistical sense. 
Lipunov et al. (1996) used for the first time a PS code 
to investigate the evolution of XRBs in the central 
part of the Milky Way over the course of 10 My (see 
Lipunov et al. 2009, for a description of their code “Sce- 
nario Machine”). A number of other codes have been de- 
veloped mostly within the last decade or so (Hurley et al. 
2002; Kiel & Hurley 2006; Belczynski et al. 2008). These 
codes have been used by several authors to carry out 
a variety of investigations, including the study of dou- 
ble compact object mergers (e.g. Belczynski et al. 2002), 
the formation of ultrashort XRBs (e.g. Belczynski et al. 
2004), the numbers and spatial distributions of XRBs 
in star clusters (e.g. Sepinsky et al. 2005), the evolution 
of XRBs in a brief star-formation episode as a model 
of starburst systems (Eracleous et al. 2006), the forma- 
tion of binary millisecond pulsars (e.g. Pfahl et al. 2003; 
Hurley et al. 2010), binary fractions in globular clusters 
(e.g. Ivanova et al. 2005; Hurley et al. 2007), and num- 
bers and birthrates of symbiotic XRBs in the Galaxy 
(e.g. Lii et al. 2012). 

Physically motivated PS modeling and detailed com- 
parisons of XLF characteristics can be used to under- 
stand how these observations are linked to the forma- 
tion and evolution of XRB populations in galaxies. This 
type of work was pioneered by Belczynski et al. (2004) 
who compared a theoretical X-ray luminosity function 
(tXLF) with observational X-ray luminosity functions 
(oXLFs) for XRBs in the dwarf irregular galaxy NGC 
1569 obtained with Chandra. Linden et al. (2009, 2010) 
studied the XLF for HMXBs and Be XRBs in the SMC. 
Fragos et al. (2008, 2009) modeled the XLFs in the two 
elliptical galaxies NGC 3379 and NGC 4278, and investi- 
gated the contributions from subpopulations of LMXBs. 
Zuo & Li (2011) used the PS code of Hurley et al. (2002) 
to investigate the X-ray-evolution of late-type galaxies 
over ~ 14 Gy of cosmic time. 

To carry out this type of work, it is necessary to com- 


bine PS models with star formation history (SFH) in- 
formation for specific galaxies. For their elliptical galax- 
ies Fragos et al. (2008, hereafter F08) assumed an initial 
J-function star formation episode. For later type galax- 
ies SFH can be obtained via spectral energy distribution 
(SED) modeling, which requires multiwavelength infor- 
mation for a given galaxy. 

In this paper we use a sample of 12 nearby galax- 
ies from the Spitzer Infrared Nearby Galaxy Survey 
(SINGS) covering a range in star forming properties to 
extend binary PS modeling to later type systems. We 
construct oXLFs for their off-nuclear point source XRB 
populations and tXLFs by combining population synthe- 
sis modeling results, obtained by means of the PS code 
StarTrack, and star formation histories from the liter- 
ature. This allows us to compare the two sets of XRB 
XLFs and investigate the range of acceptable values for 
XRB formation and evolution parameters. 

This paper is part of a larger effort to understand 
the formation and evolution of extragalactic XRBs by 
means of the most advanced PS modeling to date. Other 
papers in the series include Fragos et al. (2013, here- 
after F13) and Tremmel et al. (2013, accepted, as- 
troph/1210.7185, hereafter T13). F13 study the evo- 
lution of the global XRB population with redshift by 
using the Millennium-II simulation as initial condi- 
tions. They accurately reproduce local group HMXB 
and LMXB luminosity scaling relations with SFR and 
M* (Lehmer et al. 2010; Mineo et al. 2012), respectively. 
T13 use the same initial conditions to explain observa- 
tional XLFs for the integrated XRB emission from entire 
galaxies (Tzanavaris & Georgantopoulos 2008) and make 
predictions for higher redshifts. In this paper we apply 
the same grid of PS models, combining it with SFH in- 
formation for nearby galaxies (Noll et al. 2009). 

The structure of the paper is as follows: In Section 2 
we present the observational sample and the construction 
of oXLFs. Section 3 discusses the calculation of tXLFs. 
Section 4 presents likelihood functions for establishing 
best PS models. Results are presented and discussed 
in Section 5. Section 6 gives a summary and discusses 
future prospects. 

2. OBSERVATIONAL SAMPLE AND XLFS 

We use galaxies selected from the Spitzer Infrared 
Nearby Galaxy Survey (SINGS, Kennicutt et al. 2003). 
This survey was designed to be a diverse sample of intrin- 
sic galaxy properties with multiwavelength data ranging 
from the ultraviolet to the far infrared. As part of a 
large Chandra program (XSINGS, Jenkins et al. 2010), 
the Advanced CCD Imaging Spectrometer (ACIS) was 
used to extend the survey’s wavelength coverage to the 
X-ray regime. Details regarding the sample selection, X- 
ray observations, source detection and characterization 
will be presented in a forthcoming publication (Jenk- 
ins et al., in prep.) Briefly, basic X-ray data reduction 
was carried out using standard Chandra X-ray Center 
tools. Point source detection was performed in the soft 
(0. 3-2.0 keV), hard (2.0-10.0 keV) and full (0.3-10.0 kev) 
band with CIAO 11 WAVDETECT to construct a candi- 
date source list. The final list was produced by using 
the software acis extract (ae, Broos et al. 2010) to 

11 http://cxc.harvard.edu/ciao 
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perform aperture photometry and produce a catalog of 
point sources with associated fluxes and luminosities for 
each galaxy. 

We use a sub-sample of 12 SINGS galaxies that have 
SFHs from spectral energy distribution fitting with the 
code CIGALE 12 (Noll et al. 2009). Galaxies are also 
selected to have at least 15 detected, non-nuclear, X- 
ray point sources as a prerequisite for the production of 
meaningful XRB XLFs. Details of the galaxy sample 
are given in Table 2. Spitzer - infrared and Chandra-X- 
ray galaxy images are shown in Fig. 1. To illustrate the 
properties of our sample relative to the rest of the SINGS 
galaxies, in Fig. 2 we show both a color-magnitude dia- 
gram and a plot of star formation rate vs. stellar mass 
for the full SINGS sample, highlighting our galaxies. In 
both plots, galaxies are separated into three broad mor- 
phological categories, namely E/SO (shown in red), Sa- 
Sbc (magenta), and Sc-Im (blue). Compared to the rest 
of the SINGS galaxies, our systems have intermediate to 
high SFR and M*, and intermediate to red colors. 

At low point source luminosities incompleteness effects 
arise, compromising the construction of unbiased XLFs. 
These can be mitigated either by limiting observational 
XLFs to the luminosity range in which incompleteness is 
not significant or by performing incompleteness correc- 
tions. The second approach is preferable since it allows 
the construction of XLFs covering a wider dynamic range 
in X-ray luminosity. We use the method of Zezas et al. 
(2007, see also Zezas & Fabbiano (2002)) to create sim- 
ulated X-ray source catalogs, calculate the source de- 
tection probabilities and obtain incompleteness correc- 
tions as a function of source and background intensity 
(in counts) and off-axis angle for X-ray detected sources. 
Note that this detection probability is otherwise inde- 
pendent of other source characteristics such as detection 
band and intrinsic source spectrum. 

oXLFs are shown in Fig. 3 by the blue curves before 
(dotted) and after (solid) completeness correction. In 
practice, we find that completeness corrections are small 
and do not change the models that are associated with 
the highest-likelihood tXLF for a given galaxy. Specifi- 
cally, a quality check on our best populated galaxy NGC 
1291, for which we apply both methods, shows that our 
conclusions do not change. We note that the highest- 
likelihood tXLF is also shown in Fig. 3 (in red and dark 
grey; see Sec. 4 for details). 

3. POPULATION SYNTHESIS MODELING 
3.1. StarTrack 

The main tool we use to perform our PS simula- 
tions is StarTrack, a state-of-the-art PS code that has 
been tested and calibrated using detailed mass transfer 
star calculations and observations of binary populations. 
StarTrack has been applied to numerous interpretation 
studies of X-ray and radio pulsar binary populations, as 
well as 7-ray bursts and binary black holes in the context 
of gravitational-wave sources (see detailed description in 
Bclczynski et al. 2008, and references therein). In sum- 
mary the code incorporates all the important physical 
processes of binary evolution: 

(i) The evolution of single stars and non-interacting bi- 

12 http://cigale.oamp.fr/ 


nary components from zero-age main sequence to rem- 
nant formation is followed with the use of high-quality 
analytic formulae (Hurley et al. 2000). Various wind 
mass loss rates dependent on stellar evolutionary stage 
are incorporated and their effect on stellar evolution is 
taken into account. 

(n) Changes in all the orbital properties are tracked. 
Through numerical integration of four differential equa- 
tions, the evolution of orbital separation, eccentricity and 
component spins is tracked; these depend on tidal inter- 
actions as well as angular momentum losses associated 
with magnetic braking, gravitational radiation and stel- 
lar wind mass losses. 

(iii) All types of mass-transfer phases are calculated: 
Stable, driven by nuclear evolution or angular momen- 
tum loss, and thermally or dynamically unstable. 

(iv) Supernova explosions are treated accounting for 
mass loss and asymmetries with natal kicks to neutron 
stars and black holes at birth; systemic velocities for all 
binaries are calculated. 

(v) The calculation of mass transfer rates in binaries with 
accreting neutron stars and black holes (driven by stel- 
lar winds or Roche-lobe overflow) has been calibrated 
against detailed mass transfer sequences, and X-ray lu- 
minosities are calculated by incorporating appropriate 
band-pass corrections and spectral models. 

This paper uses results that are based on a recent ma- 
jor revision of the StarTrack code that includes updated 
stellar wind prescriptions and their re-calibrated depen- 
dence on metallicity (Belczynski et al. 2010). Two newer 
updates have not been taken into account, as our sim- 
ulations were complete before these updates had been 
implemented. For reference, these are (1) a revised neu- 
tron star and black hole mass spectrum, leading to fully 
consistent supernova simulations (Belczynski et al. 2012; 
Fryer et al. 2012); and (2) a more physical treatment 
of donor stars in common envelopes via actual A val- 
ues (Dominik et al. 2012), where A is a measure of the 
donor’s central concentration and the envelope binding 
energy. 

Table 1 lists the full set of input parameters used in our 
PS modeling. We construct a large PS model grid by us- 
ing a range of values for these parameters as indicated in 
the table. For detailed discussion of how each parameter 
affects the overall XRB population, we refer the reader 
to F13 (Section 5.1 and Fig. 6) and T13 (Section 4.3 and 
Fig. 7). 

The input parameters fall into two categories. First, 
parameters that mostly characterize initial properties of 
the population, such as initial mass function (IMF), ini- 
tial binary mass ratio (g), and distribution of initial or- 
bital separations. For these parameters we usually have 
information from observational surveys of binary stars, 
constraining the range of values used. The second cate- 
gory comprises parameters that are associated with phys- 
ical processes that are poorly understood, such as the 
efficiency, «ce, of converting orbital into thermal energy 
that will be used to expel the donor’s envelope during a 
common envelope (CE) phase. In this work the related 
parameter varied is Aocej which is the product of this 
efficiency and the donor central concentration, A. 

The complete grid of 288 PS models is the same as 
that used by F13 and T13. To construct this grid, 
we vary all parameters known from earlier studies to 
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affect XRB evolution and formation of compact ob- 
jects (Belczynski et al. 2007; Fragos et al. 2008, 2009; 
Linden et al. 2009; Belczynski et al. 2010; Fragos et al. 
2010). Specifically we use: 

• Four A«ce values (0.1, 0.2, 0.3, 0.5); 

• Three stellar wind strengths, rywind (0.25, 1.0, 2.0). 
This parameter is used to multiply the stellar wind 
prescription of Belczynski et al. (2010); 

• The distribution of natal kicks for BHs formed 
through direct collapse (no kicks or 10% of the 
Hobbs et al. (2005) distribution for NSs); 

• A CE-HG flag for systems with a donor in the 
Hertzsprung gap, either allowing all possible com- 
mon envelope events or always imposing a merger 
(Belczynski et al. 2007); 

• Three distributions of binary initial mass ratios, 
q — T/ sef:OTK ];.jry /i\ / pvirmiry ■ This distribution speci- 
fies the mass of the secondary, whereas the mass of 
the primary is governed by the IMF. For our full 
model grid, we use a uniform (flat) distribution, 
q = 0 — > 1, a twins distribution, q = 0.9 —> 1, 
and a mixed distribution, 50% uniform and 50% 
twins. However, in this paper we do not use the 96 
twin q models 97 — 192, as F13 have already clearly 
shown that these are very inadequate in reproduc- 
ing observed XRB populations, as they prevent the 
production of LMXBs altogether. This limits the 
models used in this paper to those in the ranges 
1 — 96 and 193 — 288, i.e. 192 models in total. 

In addition, for each of the 288 models in our full grid, 
we also use nine metallicities, keeping all other param- 
eters same. In this paper we only consider solar metal- 
licity models, as the best estimates of Moustakas et al. 
(2010) using two different methods appear to straddle so- 
lar metallicity for all SINGS galaxies in our study. The 
best fit SEDs that we convolve with StarTrack models 
to construct tXLFs for individual galaxies also assume 
solar metallicities (Noll et al. 2009). 

Each model follows the evolution of 5.12 x 10 6 stars 
over 14 Gyr. Since we are only using models 1 — 96, 
corresponding to a uniform q distribution, and models 
193 — 288, corresponding to a 50% - 50% mixed q distri- 
bution, in relevant figures we indicate these model ranges 
for clarity. For reference, the full list of 288 models and 
associated parameters can be found in Table 4 of F13. 

3.2. Theoretical X-ray luminosities 

Since our goal is to construct XRB XLFs, we need 
to identify all binaries in our simulations that become 
XRBs and register their X-ray luminosities as a func- 
tion of time. XRBs are mass-transferring binary stel- 
lar systems with a compact object accretor, either a 
black hole (BH) or a neutron star (NS). Systems with 
donors less massive than 3 M & are labeled LMXBs, and 
vice versa for HMXBs. LMXBs are always Roche lobe 
overflowing (RLOF) systems, while HMXBs are usually 
wind-fed, but can also exhibit RLOF behavior. Accord- 
ing to whether they undergo thermal disk instability or 


not, RLOF systems can be either transient or persistent, 
while wind-fed systems are always persistent. 13 
We follow the methodology described in F08 and 
Fragos et al. (2009) to identify all model XRB sources 
and keep track of properties essential for estimating their 
X-ray luminosity, Lx, as a function of time. These prop- 
erties include the mass- transfer rate, M , as well as the 
mass and radius of the accretor, M a and R a , respec- 
tively. We also identify BH or NS accretors, transient 
or persistent sources, evolutionary stages of donors and 
donor masses. 


3.2.1. RLOF systems 

For persistent systems, Lx is estimated as 


L x = min ^/L E dd, ??boie 


GM a M 

Ra 


( 1 ) 


with / equal to unity. The value of R a is 10 km for a 
NS and 3 Schwarzschild radii for a BH, e gives a conver- 
sion efficiency of gravitational binding energy to radia- 
tion associated with accretion onto a NS (surface accre- 
tion, e = 1.0) or BH (disk accretion, e = 0.5), and ?7boi 
is a factor that converts the bolometric luminosity to 
Lx,o.3-io.okeV, the X-ray luminosity in the full Chandra 
energy band, consistent with our observations. We use 
results from the literature to tabulate the best available 
estimates of this factor and its la uncertainty for neutron 
star and black hole accretors in “low-hard” and “high- 
soft” state systems (see F13 for a detailed discussion and 
references). Since we are interested in combining our PS 
models with galaxy star formation histories, we define 
lookback-time windows, St ib , in which we evaluate the 
total stellar mass produced in the galaxies (see Sec. 3.3). 
In the StarTrack simulations each model source is asso- 
ciated with its own time-step window, St s t ep , which we 
use to calculate the fractional time, 6N, that a source is 
on in a given galaxy lookback time interval, i.e. 


SN = 6t step /5t lh . (2) 

For transient systems we follow the prescription of F08, 
calculating the outburst X-ray luminosity depending on 
the accretor type, either BH or NS, via 

( min (/bh x L E dd,/BH x ^Edd [i^])jf° r BH 

Lx = ^bolCX < 

{ min (/ns x i E dd, G fl ° A /° rit ) 1 f° r NS 

(3) 

Here /bh, /ns are factors setting an upper limit for the 
maximum Eddington luminosity for black holes and neu- 
tron stars, equal to 2 and 1, respectively. Merit, -Mrf and 
P are the critical mass transfer rate for thermal disk in- 
stability, the rate at which the donor star is losing mass 
and the orbital period, respectively. In addition, tran- 
sient systems are mostly in a quiescent state and are too 
faint to be detectable, except when they go into out- 
bursts. The fraction of time they are in outburst defines 

13 Be XRBs are a special case, since, although they are wind- 
fed, they show quasi-periodic outbursts. However, the origin of 
this behavior is not the thermal disk instability. 
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their duty cycle, DC. Following Fragos et al. (2008), we 
estimate DC as (Dobrotka et al. 2006) 


DC = 



( 4 ) 


for NSs and 5% for BHs (Tanaka & Shibazaki 1996). For 
transient systems it then follows that the fractional time 
a source is on in a galaxy lookback time interval becomes 
5N = Ststep/dtib x DC. 


3.2.2. Wind-fed. systems 

For wind systems the Lx calculation is as for for per- 
sistent RLOF systems (Equation 1), while the fractional 
lookback time is once more given by Equation 2. 

3.3. Stellar Mass Normalization 

The StarTrack simulations described thus far link stel- 
lar mass described by a set of physical parameters with 
an XRB population but do not contain information to 
link the XRB population to host galaxy properties. Thus 
the stellar mass produced by a given model is in a sense 
arbitrary. To construct tXLFs corresponding to real 
galaxies, we need to modify this arbitrary mass by mak- 
ing use of the galaxy’s star formation history. This in 
turn modifies the numbers of XRB sources produced by 
the model. For the galaxies in our sample, SFH estimates 
exist based on SED fitting by Noll et al. (2009, see Ta- 
ble 3). We use these results to calculate the total stellar 
mass produced for each SINGS galaxy in our sample in 
each lookback time window. Noll et al. have used the 
SED fitting code CIGALE and assummed two exponen- 
tially varying star formation rates for a young and an 
old population. Using their symbols (Table 3), the total 
SFR at a lookback time fib is given by 


SFR tlb = 


/burst Afg a l 


TySP 


otySp/fySP _ 


1 ) 


D *lb/PySP 


(1 /burst )Afgal c t lh / ToS 

ToSP (e toSp / T ° SP - l) 


( 5 ) 


As in Table 3, M ga i is the total mass of stars and gas 
that originates from stellar mass loss in the galaxy, /burst 
is the mass fraction of the young population, r y sp, r Q sp 
are the e-folding timescales for the young and old stellar 
population, and t y sp , f 0 SP the ages of the young and old 
stellar population. All of these parameters are SED fit 
results allowing us to calculate SFRt lb via equation 5, 
and then stellar masses, used to scale XRB numbers for 
each galaxy-model pair. 


3.4. Construction of theoretical XLFs for SINGS 
galaxies 

As explained in F13, the bolometric corrections used to 
convert StarTrack-derived Lx values for model sources 
to the full Chandra band are empirical and introduce 
an uncertainty to the estimated X-ray luminosity value. 
Thus each Tx,o.3-io.okeV value can be thought of as the 


mean of a Gaussian distribution with a standard devi- 
ation originating in the uncertainty introduced by the 
bolometric correction. 

Further, as explained above, for each model source the 
fractional time a source is on in a galaxy lookback time 
window is given by SN. This number can also be consid- 
ered as the mean of a Poisson distribution that represents 
the expected number of times that this source will ap- 
pear in the XLF. For each source, we use this information 
to draw a random Poisson deviate, giving a number of 
times that this source will be on. For each case that the 
source is on, we draw a random Gaussian deviate from 
the Lx distribution giving an Lx value for that appear- 
ance. When this procedure is complete for all sources, we 
obtain a set of Lx values, which can be converted to an 
XLF. However, this is a single realization. To obtain a 
reliable estimate for the mean XLF and its uncertainty, 
we carry out 500 Monte Carlo realizations of this pro- 
cess for the total XLF, and 100 for subpopulation XLFs 
(donor/accretor type, LMXB/HMXB). This procedure 
gives a reliable estimate for the mean number of sources 
in each Lx bin, as well as ±ler uncertainties. 

Overall, we thus obtain 12 x 192 = 2304 tXLFs (in 
both cumulative and differential form) combining SFHs 
for 12 SINGS galaxies and 192 StarTrack models. The 
tXLFs for the best models (Sec. 4) are shown as red 
dotted curves in Fig. 3. 

3.5. Background contamination 

Our purpose is to compare our theoretically derived 
tXLFs for XRB populations to oXLFs of point sources 
in SINGS galaxies. Whereas all theoretically obtained 
point sources in the tXLFs are known to be XRBs by 
construction, this is not necessarily the case with all ob- 
served point sources in SINGS galaxies. It is possible 
that some of the latter are actually background AGN 
rather than galactic XRBs. We correct our pure XRB- 
based tXLFs by adding a component that takes into ac- 
count the expected number of background X-ray point 
sources in the area surveyed for each galaxy. 

We correct our differential and cumulative theoretical 
XLFs as follows. We use the log IV — log S' results of 
Kim et al. (2007) (r = 1.7, broad band B, their Table 
3) to estimate the number of background sources ex- 
pected in each luminosity bin of our theoretical XLFs, 
taking into account the area observed for each galaxy 
(either the Chandra S3/I0-I4 detector area or the D 25 
region, whichever is smaller). We then add these num- 
bers to our “pure XRB” source numbers in each lumi- 
nosity bin, thus obtaining background corrected tXLFs. 
In Fig. 3 these corrected, final tXLFs are shown as solid 
red curves for highest-likelihood models (Sec. 4). A com- 
parison with observational completeness corrections (see 
the blue curves in Fig. 3), which only affect the faint end 
of the XLF, shows that the background correction affects 
the XLF over the full range of X-ray luminosities. 

4. LIKELIHOOD FUNCTIONS 

To obtain a quantitative estimate of the level of agree- 
ment between observational and theoretical XLFs, we 
construct and evaluate likelihood functions, using the dif- 
ferential XLF versions. As explained below, we calculate 
two types of likelihoods. 
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4.1. Pair likelihoods for a given galaxy 

For each galaxy there is a unique, completeness cor- 
rected oXLF coming from our Chandra observations 
(continuous blue curves in Fig. 3). We wish to establish 
which of the 192 StarTrack models used in this work, af- 
ter it has been combined with a specific galaxy’s SFH and 
corrected for background AGN, best describes this oXLF. 
In other words, there are 192 tXLFs for this galaxy that 
need to be compared with a single oXLF. Thus, we cal- 
culate pair likelihoods for 192 o-t XLF pairs. 

Given a galaxy and its oXLF, k, and a model tXLF, 
m, we define the k — m th pair likelihood 1 £ pa ir,fcm, as the 
probability, P, of obtaining an observational set of XLF 
data points, k, given a theoretical model XLF, m. This 
is given by 

^pair ,km — n Ppoisson(A^ 0 bs,i; (6) 

i 

Here Pp 0 i SS on(N ohSii , N thyi ) is the Poisson probability 
of observing iV 0 bs,i point sources, in the i th luminosity 
bin, treating the theoretically obtained number of point 
sources, N t h,i, in the same bin, as the expectation value 
for the number of sources at this luminosity. The total 
pair likelihood is the product of all such probabilities over 
all luminosity bins. This compares the agreement of the 
two XLFs in each luminosity bin, and thus its overall 
shape and normalization. 

Note that tXLFs are calculated for 100 equal-sized bins 
in log Lx , spanning the observational range of luminosity 
values. The oXLFs are binned to match the tXLFs bins. 
In addition, in order to compare corresponding quanti- 
ties, ./Vobs.i are the observed numbers before correction for 
incompleteness, while N t h,i are the theoretical numbers 
after the observationally derived incompleteness informa- 
tion for a given bin has been taken into account. 

For a given galaxy, the best model is that for which this 
procedure produces the highest pair likelihood value. We 
tabulate highest-likclihood results based on this proce- 
dure in Table 4 for all galaxies in our sample. Since the 
absolute numeric value of £ pa ir,fcm obtained via Equa- 
tion 4.1 has no particular meaning, in column 10 each 
£ pa ir ,km value is normalized to the highest value in the 
table. This leads to a ranking, shown in column 9, with 
the highest £ pa i r ,fcm value having rank 1. We also plot 
likelihood results against model number for all o-t pairs 
(i.e. not just highest-likelihood pairs) in the top panel of 
Fig. 4. For convenience, in this figure £ pa i r ,fcm values are 
normalized so that they range from 0 to 1 (see Sec. 4.3). 

4.2. Global likelihoods 

The previous procedure determines the best model for 
the k th oXLF by estimating the probability, £ pa ir ,km, of 
obtaining this XLF with each model. We are also inter- 
ested in knowing the best model, m, for all oXLFs taken 
as a whole. To determine this, we calculate the product 
of all £ pa ir,fcm values by defining the global likelihood 

-^global, m = | | T pa j r ./rm (7) 

k 

where £ pa ir,fcm is given by Equation 4.1 and the index k 
runs from 1 to 12, corresponding to each of the 12 oXLFs. 
The results of this procedure are tabulated in Table 6, 


ranked by the ratio of each global likelihood value to 
the maximum likelihood value in the table (model 245). 
In Fig. 3 we also show for each galaxy the tXLF that 
corresponds to the best global model (245) as a dark 
grey curve. 

4.3. Likelihood Normalization 

The specific pair or global numeric likelihood values 
obtained have no particular meaning, except relative to 
each other. For plotting purposes (Figures 4, 5) it is con- 
venient to normalize likelihood values so that the max- 
imum value is 1 and the minimum 0. This is done as 
follows. 

Let C stand for either £ pa j r or £ g i 0 b a i- Then the final 
normalized likelihood value is given by 

An = A m i n / max( A m ; n ) (8) 

where 

A min = ln[£/min(£)] . (9) 

In what follows we use an additional subscipt to specify 
whether the final normalized likelihood is a pair or a 
global likelihood (A niPa j r and A n>g i 0 b a i, respectively). 

5. RESULTS AND DISCUSSION 
5.1. Best models 

A visual inspection of o-t XLF pairs in Fig. 3 sug- 
gests that in many cases tXLFs (solid red curves) are 
successfully reproducing oXLFs (solid blue curves) both 
in shape and normalization. However, there is consider- 
able variation and room for improvement. The best case 
is NGC 4826 and Model 269, which based on its likeli- 
hood estimate has rank 1 in Table 4. The tXLF for this 
galaxy and Model 245 (our best global model) , shown by 
the dark grey line in the panel, is also very close to the 
oXLF. The worst case is NGC 3184 and Model 277 (rank 
12 in Table 4). Since likelihood estimates take into ac- 
count source numbers in luminosity bins (Equation 4.1), 
the low likelihood estimate in this case is driven by the 
strong discrepancy in numbers at low luminosities. In 
many cases XLFs for the best pair model (red curve) 
and best global model (dark grey curve) are very close, 
but there are also cases where these are in strong dis- 
agreement. 

The top panel of Fig. 4 displays the normalized pair 
likelihood, A n pa i r , for all 192 x 12 o-t pairs (all 12 galax- 
ies), versus all models used in this paper. The imposed 
normalization (see Sec. 4.3) is such that the maximum 
normalized likelihood for each o-t pair is 1 and the min- 
imum is 0. Curves plotting normalized likelihood values 
as a function of model number (1-96 and 193-288) for 
each galaxy/o-t pair are shown with a different color. 
It is immediately obvious that for each galaxy there are 
a number of good models, indicated by maxima in the 
curves, as well as a number of poor models indicated 
by minima. Closer inspection of the plot reveals the re- 
markable fact that, overall, there is strong clustering of 
local maxima and minima, indicating that good models 
are good for all galaxies (and vice versa for poor models). 
This is so despite the fact that for a given galaxy the best 
models do not always match oXLFs well (Fig. 3). This 
suggests that, to first order , the parameters associated 
with a good model are useful indicators of XRB physics 
across all 12 galaxies. 
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This global trend is corroborated by the lower panel 
of Fig. 4 which shows the normalized global likelihood, 
A n , global- As above, the normalization is such that the 
maximum is 1 and the minimum 0. Given the similarity 
of A n pa i r curves (top panel) with each other, it is not 
surprising that they are also similar with the global like- 
lihood curve, A n global: as the latter combines, for a given 
model, all pair likelihoods for all galaxies (Equation 7). 

To understand what highest-likelihood results imply 
for individual model parameters, we tabulate results in 
two ways. We first rank highest-likelihood galaxy-model 
pairs according to pair likelihood value and show the re- 
sults in Table 4. Clearly, this tabulation favors a Actce 
value of 0.1, a high-end IMF exponent of —2.7 and a 
mixed initial q distribution. Second, we rank all 192 
models according to global likelihood value and show the 
results for the 15 highest ranked models in Table 6 14 . 
These correspond to the 15 highest peaks in the lower 
panel of Fig. 4. These results also favor a A«ce value of 
0.1, and to a lesser extent a high-end IMF exponent of 
—2.7 and a mixed initial q distribution. 

By comparing Tables 4 (col. 2) and 6 (col. 1), we note 
that all six best global models are also among the best 
pair models. This just underlies the fact that these are 
the best models overall. Conversely, nine of the best pair 
models are also among the best global models. For the 
three galaxies NGC 1291, NGC 2841 and NGC 4736, the 
best pair models are not among the 15 best global mod- 
els. It is, however, remarkable that for nine out of twelve 
galaxies one of the best 15 models that represent global 
averages over all galaxies also describes the individual 
galaxy XLF. 

We note that for galaxy NGC 1291 Luo et al. (2012) 
use Chandra observations that are deeper than ours to 
perform a detailed study of the point source popula- 
tion. They show that high-luminosity (> 5 x 10 38 erg s _1 ) 
LMXBs in NGC 1291 are likely associated with a younger 
stellar population in the galaxy’s ring. Their deeper ob- 
servations allow them to separate the bulge from the ring 
population. Even so, we note that our results still sup- 
port their conclusions, as we find strong LMXB contri- 
butions both from our old and young populations. 

5.2. Deviations between models and observations 

A key reason for deviations of oXLFs from tXLFs is 
the presence of two “jumps” at two different luminosi- 
ties (~ 10 37 erg s -1 and a few xlO 38 erg s _1 ), which 
are not expected from observations. There are two rea- 
sons for this. First, in our StarTrack models we identify 
transient and persistent sources by comparing the calcu- 
lated mass-transfer rate of the binary to a critical mass- 
transfer rate, below which the thermal instability devel- 
ops, giving rise to transient behavior. Furthermore, we 
strictly limit the accretion rate to the Eddington limit. 
These limiting mass-transfer rates in our modeling are 
responsible for the jumps in tXLFs. However, in nature 
there is no sharp transition between thermally stable and 
unstable disks, nor a precise limit to the highest accre- 
tion rate possible. Accretion onto a compact object is 
a non-linear and much more complex process, which in 
reality can result in a smoother luminosity distribution 
with no sharp transitions. Second, the assumption of so- 

14 The full table for all 192 models is available online. 


lar metallicity in practice imposes a maximum BH mass 
of ~ 15 Mg. If for instance part of the stellar pop- 
ulation had a lower metallicity (e.g. 30% solar) then 
the maximum BH mass would increase to 30 Mq or 
more (Belczynski et al. 2010), which could smooth out 
the jump at the very luminous end of the XLF. 

More generally, it must be stressed that the SFHs we 
use are very simple. Noll et al. (2009) note that, as SED 
fitting is computationally intensive, they select a limited 
set of values for their SED parameter model grid. Thus 
the age of the old population is a constant (10 Gy), and 
there are only two possible values for the age of the young 
population, 200 Myr and 50 Myr. The parameter that 
was allowed to vary the most is “/burst”, the mass frac- 
tion of the young stellar population at the present time 
(nine possible values). Further, their detailed analysis 
of their SED fitting results shows that best-fit e-folding 
times and the age of the young population are highly 
uncertain, due in turn to photometric errors and uncer- 
tainties in stellar population models. Finally, as already 
mentioned, their SFHs have uniformly solar metallicities. 
Although this is based on the best estimates to date, 
note that these are based on gas metallicities with their 
own set of significant uncertainties (see Moustakas et al. 
2010, where two sets of metallicity results are presented, 
straddling solar metallicity for all galaxies in this sam- 
ple). Further progress in matching oXLFs will require 
more detailed SED fitting. 

5.3. Constraining Parameters 

Table 4 and especially Table 6 suggest likely best values 
for model parameters in our model grid. As mentioned, 
both tables suggest Actce — 0.1. In Fig. 5 we inves- 
tigate this further by plotting resistant mean 15 values 
for the normalized global value, A nig i 0 b a i, vs. the four 
StarTrack AacE parameter values used in our simula- 
tions. The calculation of the resistant mean iteratively 
rejects outliers beyond 3cr. This is useful for highlighting 
the clustering of likelihood values, which often show a 
considerable spread. 

Fig. 5 shows that higher likelihood values are sys- 
tematically favored for models with ActeE — 0.1: Not 
only do A njg i 0 b a i values converge to their highest value 
as AacE —> 0.1, but they also do so with a decreasing 
spread, as indicated by the 3cr resistant mean error bars. 

Table 4 and Table 6 also appear to favor a value of 
—2.7 for the high-end slope of the IMF. Although a value 
of —2.7 agrees with some estimates (e.g. Scalo 1986), 
we note that it is somewhat steep compared to the so- 
called “canonical” value of —2.3 established observation- 
ally for resolved stellar populations in the local group 
(Bastian et al. 2010; Kroupa 2012). In spite of the fact 
that (1) variations with metallicity and star formation 
rate density for the integrated galaxy IMF remain hotly 
debated (Kroupa 2012, and references therein) and (2) 
taking uncertainties in the canonical value into account, 
a value of —2.7 is still close to the canonical upper limit, 
we do not consider a value of —2.7 to be a robust result. 
This is because, first, the Noll et al. SED fitting results 
are assuming the canonical value, and it is these results 
that have been convolved with our StarTrack models. 

15 http:/ /idlastro. gsfc.nasa.gov/ftp/pro/robust/ 
resistant_mean.pro 


Tzanavaris et al. 


Second, four out of the fifteen best global models do, in 
fact, agree with the canonical value. 

For the remaining parameters these tables show that 
the results are inconclusive. In addition, for most of 
these parameters we only use two different values in 
the StarTrack simulations, so we cannot investigate any 
trends such as those for Aoce in Fig. 5 . 

Given the caveats for the IMF slope, the results for 
Aoce — 0.1, and to some extent for the prevalence of a 
mixed initial q distribution provide the most significant 
constraints for binary star parameters from this work. 

As A«ce is a combination of two parameters, we are 
unable to set any constraints on either A or oce indi- 
vidually. In their investigation of Galactic merger rates 
for compact objects, Dominik et al. (2012) set acE = 1-0 
and study the behavior of A. The latter is not consid- 
ered constant throughout the evolution of the donor, but 
depends on donor parameters such as mass, radius and 
evolutionary stage. They find A values between ^0.1 and 
~ 0.2 for NS progenitors, and below ~ 0.1 for BH progen- 
itors. With the assumption of acE = 1.0 our result for 
Aoce would thus be largely consistent with the detailed 
stellar evolution models of Dominik et al. (2012). 

5.4. Comparison with F13 and T13 

F13 and T13 use the same StarTrack models as 
this paper but combine them with galaxy information 
from the Millennium II simulation and the semi-analytic 
galaxy catalog of Guo et al. (2011). While F13 investi- 
gate the total galaxy specific X-ray luminosity (Lx/SFR, 
Lx/M*) evolution out to 2 ~ 20, T13 construct galaxy 
XLFs for comparison with observations out to z ~ 1.4 
and make predictions out to z ~ 20. Thus our three pa- 
pers examine XRB formation and evolution in different 
contexts, and also use different likelihood formulations. 
To investigate whether there is agreement in highest- 
likelihood models between the three papers, in Table 6 
we also show ranks from F13 and T13 for our best 15 
global likelihood models. There is good agreement over- 
all, and in certain cases the agreement is exceptionally 
good. In particular, our best ranked model (245) is also 
F13’s reference model (their rank 1) and has rank 4 in 
T13. Models 245, 277, 229, 205, 269, 249, 273, 37 and 85 
are all within the 15 highest ranked in all three papers. 
This is further evidence in support of Actce —0.1 both 
in the local universe and over cosmic time. 

5.5. XLF subpopulations 

In Figures 6 and 7 we plot the total tXLF for each 
galaxy and highest-likelihood model (black continuous 
curve), together with constituent subpopulation tXLFs. 

We define two sets of subpopulations. The first is 
shown in Fig. 6 where we plot XLFs for LMXBs (red) and 
HMXBs (blue). In our simulations XRBs are labelled 
HMXBs if the donor star has mass M<j on0 r > 3 Mq, and 
LMXBs otherwise. The mixture of populations present 
in each galaxy is closely tied to the assumed SFH (Ta- 
ble 3). To illustrate this, XRBs originating in the “old” 
stellar population are shown by the dashed curve in 
Fig. 6, while those originating in the “young” SP are 
shown by the solid curve. Note however that the terms 
old and young can be misleading in this context. Al- 
though the two are distinct in terms of age, with the 


young one appearing several Gyr after the old one, it is 
possible for the latter to still be actively star-forming, de- 
pending on the e- folding timescale r Q gp. In other words, 
an “old” population is not necessarily a “red and dead” 
one. 

The trends seen in Fig. 6 for each galaxy can be un- 
derstood qualitatively by referring to (1) the SFH (Ta- 
ble 3), which we combine with StarTrack models to con- 
struct tXLFs, and (2) Table 4 which gives the details of 
each best model for each galaxy. We should also keep 
in mind the relative evolutionary timescales for LMXBs 
and HMXBs. As shown in F13 16 (their fig. 2), in a single 
burst population the contribution in X-ray output from 
HMXBs peaks at about 5 Myr, remaining important up 
to an age of ~ 100 — 300 Myr, while LMXBs take over 
at ~ 100 — 200 Myr. In terms of the assumed “young” 
and “old” SPs, a population’s contribution will be more 
significant as its e-folding timescale is longer and age is 
younger. A higher young mass fraction will tend to in- 
crease the contribution from the young population. 

Looking at galaxies NGC 1291 and 2841, we note that 
they have all best fit SFH parameters the same, and 
slightly different M ga i, while the associated StarTrack 
model is the same. It is evident that the subpopulation 
XLFs are similar, with the old LMXB SP dominating in 
both cases. 

For galaxies NGC 3184 and 3627, the main SFH dif- 
ference is a higher t q sp for the latter. Based on this, 
one might expect that in NCG 3627, the old population 
would be more dominant. However, young SP HMXBs 
dominate for most of the XLF. This is due to StarTrack 
model 245 (NGC 3627) having half the stellar wind 
strength compared to model 277 (NGC 3184). Weaker 
stellar winds lead to smaller mass loss for primaries that 
eventually become compact objects, and thus to numer- 
ous and more massive BH XRBs. The latter in turn 
tend to be more luminous than NS XRBs as (1) they can 
form stable RLO XRBs with massive companions, and 
(2) they show higher accretion rates due to the high BH 
masses. Although a weaker stellar wind also decreases 
the accretion rate in wind-fed HMXBs, it turns out that 
this is not the dominant effect (F13). 

NGC 3198 has by far the highest t 0 sp among all galax- 
ies and this is reflected in the dominance of the old 
LMXB population. In contrast, NGC 3521 and 5055 
have a much lower r 0 gp (with other parameters except 
for M ga i identical for all three) and are dominated by 
young SPs. NGC 4631 is also similar, with a somewhat 
lesser contribution from the old SP (higher r y sp). Com- 
pared to NGC 4631, NGC 4826 has lower r y sp and /burst 
leading to a relatively stronger contribution from the old 
SP. 

NGC 4736’s young SP has the smallest best-fit age (50 
My) together with a high T y gp. This is consistent with 
the clear dominance of young HMXBs, which are closest 
to their peak activity timescale (Shtykovskiy & Gilfanov 
2007), while LMXBs haven’t yet had time to contribute 
significantly to the total X-ray luminosity. 

Finally, NGC 5474 has higher t 0 sp but a higher young 
mass fraction compared to, e.g. NGC 3521, so it is mod- 
erately dominated by the young SP. 

16 F13 use the same definition for LMXBs and HMXBs as this 
paper in terms of a threshold donor mass of 3 Mq . 
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The second subpopulation set is shown in Fig. 7 in 
terms of donor and accretor stellar type. We classify 
donors as MS, evolved or degenerate. Evolved types 
have left the MS and can be in any of a number of differ- 
ent stellar evolutionary stages, including the Hertzsprung 
gap, the red giant branch (GB), core helium burning, 
early asymptotic GB, thermally pulsating asymptotic 
GB, helium MS, helium HG or the helium GB. Degen- 
erate types are helium or carbon/oxygen white dwarfs. 
On the other hand, accretors are either BH or NS. 

A number of trends are visible in these plots. Consid- 
ering the accretor subpopulations, indicated by different 
linestyles in Fig. 7, we notice that in general BH accre- 
tors (solid lines) dominate the XLFs at high luminosities. 
In addition, galaxies with no sources at Lx > 3 x 10 38 
have no BH accretors. These galaxies are NGC 3184, 
3351, 3521, 4631 and 5055. The reason for this is that, 
since BH accretors originate in the high-mass end of the 
IMF, from a statistical point of view on average there will 
be fewer such systems. In addition, BH XRBs can have 
higher luminosities than NS XRBs, but not vice versa. 
Thus, BH XRBs can only be registered only when very 
luminous sources are present (e.g. see Luo et al. 2013, 
for NGC 4649), as NS XRBs will always dominate at 
low luminosities. As a result, for galaxies with very few 
or no sources above ~ 10 38 erg s -1 models are unable to 
constrain the BH XRB population. 

All donor subpopulation XLFs, indicated by different 
colors in Fig. 7, are dominated by MS donors at low lu- 
minosities (< 10 37 erg s -1 ) 17 and by evolved donors at 
high luminosities. This is a consequence of the fact that 
evolved donors are mainly giants in RLOF systems, so 
they all have large accretion disks and long orbital pe- 
riods. Since Merit depends on the size of the accretion 
disk and orbital period, it follows that in the case of these 
systems it will always be higher than a threshold value, 
which corresponds to Lx ~ 10 37 erg s _1 (King et al. 
1996; Dubus et al. 1999).. Further, systems with mass 
transfer rates less than M cr ; t and, thus, with luminosi- 
ties lower than ~ 10 37 erg s _1 will be transients, which 
are mostly quiescent. Thus in practice at Lx > 10 37 the 
dominant donor systems will be persistent systems with 
evolved donors. 

5.6. Comparison with Lehmer et al. 2010 

It is well established that emission from HMXBs and 
LMXBs correlates with galaxy- wide SFR and M*, re- 
spectively. This is due to the fact that the former are 
relatively young (< 100 Myr) compared to the latter 
(> 1 Gyr). Lehmer et al. (2010, L10) parametrize the 
contribution of LMXBs and HMXBs to the 2 — 10 keV 
integrated luminosity in star-forming galaxies as 

Lhx = £hx XB + ^hx XB — a + /3 SFR , (10) 

where Lhx = Lx, 2 .o-io.okeV- They use a sample of 
galaxies spanning a large range in star-forming activ- 
ity observed with Chandra to obtain best-fit values a = 
(9.05 ± 0.37) x 10 28 erg s^Mq 1 and j3 = (1.62 ± 0.22) x 
10 39 erg s -1 (MQyr -1 ) -1 . We use these results to com- 
pare with our work, noting also that they are consis- 
tent with other work both on the Lx—M„ relation in el- 

17 This is not always obvious in Fig. 7 as the Lx ranges shown 
are aimed to match the observed XLF ranges. 


liptical galaxies (Gilfanov et al. 2004; Zhang et al. 2011; 
Boroson et al. 2011) and on the Lx~ SFR relation in 
galaxies with high specific SFR, where HMXBs are dom- 
inant (Gilfanov 2004; Mineo et al. 2012). 

We use our subpopulation XLFs for LMXBs and 
HMXBs to compare with the L10 results as follows. 
We obtain galaxy-wide X-ray luminosities by integrat- 
ing each model XLF over all luminosity bins. Note that 
for each galaxy we use both the highest pair-likelihood 
model XLF and the XLF which is based on the best 
global model 245, thus obtaining two sets of galaxy-wide 
X-ray luminosities. We then use the M* and SFR val- 
ues for each galaxy to estimate Lhx XB an d T B x XB from 
Equation 10. Since our data are in the 0.3 — 10 keV 
band, we convert Lrx values to Lx,o.3-io.okeV by using 
the mean Lx,o.3-io.okeV to Lhx ratio based on the F13 
models. In Fig. 8 we plot the integrated X-ray luminosi- 
ties for LMXBs and HMXBs against M* and SFR. The 
grey open circles are for integrated luminosities that use 
the highest pair-likelihood model XLF, while the filled 
black circles are for model XLFs using the best global 
model 245. We also show the relations M* and 

Thx XB ~ SFR for the best-fit a, /3 values of L10, with 
the associated ler scatter. 

The comparison shows that, when the highest pair- 
likelihood models are used, there is some agreement, es- 
pecially for LMXBs, but some of these these fail to repro- 
duce the L10 relations. In contrast, when the best global 
model is used, integrated luminosities are in much better 
agreement with the L10 relations, at least for HMXBs. 
The residual scatter is likely due to limitations of SFHs 
and uniformly solar metallicities. The observed agree- 
ment suggests that averaging over all twelve galaxies to 
calculate global likelihoods produces results that, in a 
statistical sense, are more reliable. On the other hand, 
individual pair-likelihoods inevitably suffer from all SFH 
uncertainties discussed earlier, even though a given tXLF 
might match a given galaxy oXLF better than the tXLF 
for model 245. This result is fully consistent with the fact 
that model 245 is also the top-ranked model of F13. Us- 
ing an independent likelihood formulation, these authors 
constrain their best models by comparing with observa- 
tional work, including the L10 relations. 

This result provides strong motivation for future work. 
One would expect that for a larger dataset, such as all 75 
SINGS galaxies, agreement would be further improved. 

6. CONCLUSIONS 

We have constructed theoretical XRB XLFs, for 
the first time corrected for background contamination 
(Sec. 3.5) and including ler uncertainties, for 12 nearby, 
late-type SINGS galaxies. We compare them to obser- 
vational XRB XLFs, corrected for incompleteness, by 
means of a likelihood approach (Sec. 4). 

Our main results are as follows: 

1. By comparing 192 theoretical models (Sec. 3.1) to 
observed XRB populations in twelve nearby galax- 
ies (Sec. 2), we are able to constrain best values 
for XRB formation and evolution parameters (Ta- 
bles 4 and 6). This is the largest scale comparison 
in terms of numbers of nearby galaxies and theo- 
retical models to date. 
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2. There is substantial range in the level of agreement 
between observational and theoretical XLFs for in- 
dividual galaxies, due to SFH and some model lim- 
itations (see Sec. 5.2 and below). For about half 
of the galaxies the agreement is not good. How- 
ever, for any given model, likelihoods are consis- 
tently high or low both when estimated for indi- 
vidual galaxies (£ pa i r ) and when averaged over the 
full galaxy dataset (£ g i 0 bai)- Thus parameters as- 
sociated with highest-likelihood models provide in- 
sight for XRB physics irrespective of the details of 
specific galaxies (Figures 4 and 5). 

3. Our best models have AacE — 0.1 and a mixed 
initial q distribution. Their IMFs have high-end 
slopes of —2.35 or —2.7 and rewind — 1-0 — 2.0 (Ta- 
bles 4 and 6, Fig. 5). However, we stress that fur- 
ther work is required before reliable values for these 
parameters can be established (see below). 

4. Our best models for XRBs in nearby galaxies are 
in agreement with work describing the cosmological 
evolution of XRBs (F13) as well as integrated XRB 
emission from entire galaxies (T13, see Sec. 5.4 and 
Table 6). 

5. Model XLFs show considerable variation in their 
constituent systems. Some galaxies have no BH 
accretors and most have a substantial contribution 
from LMXBs (Sec. 5.5). 

6. The integrated model XLF X-ray luminosity due 
to LMXBs and HMXBs that is based on the best 
global model in this paper and in F13 agrees with 
the expectations from the L10 relations based on 
galaxy-wide stellar masses and, especially, SFRs 
(Sec. 5.6 and Fig. 8). 

This paper represents the first conserted effort to 
model observational XRB XLFs for a set of late- type 
galaxies of this size. For individual galaxies, tXLFs 
match oXLFs with varying degrees of success. Even 
so, there are clear global trends regardless of individ- 
ual galaxies. We can thus begin drawing conclusions for 
a number of physical parameters related to XRB forma- 
tion and evolution. 

The major limitations for this work come from the 
SFH, the small sample, imposed limiting mass-transfer 
rates (Sec. 5.2) and poor understanding of many aspects 
of physics related to XRBs, precluding the construction 
of good models. Observationally, this provides motiva- 
tion for increasing the sample to include more nearby 
galaxies with reliable SFHs. On the computational side, 
work on the physics of XRB formation and evolution 
needs to include detailed modeling of the common enve- 
lope phase, as well as detailed self-consistent mass trans- 
fer calculations. Poor or non-existent understanding of 
physics remains a challenge for models. Thus, we do 
not really know how either BH or Be XRBs form, al- 
though the latter constitute an important population 
which forms the majority in the Small and Large Magel- 
lanic Clouds. We also do not yet understand the physics 
of disk instability. 

Chandra 1 s superb angular resolution is critical for this 
type of work. Deeper Chandra observations would al- 
low to expand the dynamic range for comparisons with 


models to fainter luminosities, mitigating the need for 
completeness corrections. 

Deep HST observations to securely identify counter- 
parts can also further our understanding of the nature of 
XRB systems. Such a task has been successfully achieved 
to date only for a handful of nearby galaxies (at dis- 
tances < 10 Mpc), uncrowded regions and the bright- 
est stars (e.g. Dalcanton et al. 2009; Rejkuba et al. 2009; 
Tikhonov & Galazutdinova 2012). In the medium term, 
the advent of 30-meter class telescopes such as the ELT 
and TMT, promises to pave the way to major break- 
throughs in this held (Greggio et al. 2012). 
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Figure 1. Infrared morphologies and X-ray point sources of SINGS galaxies in this paper. Each thumbnail pair shows a Spitzer - IRAC 
composite false-color image at left (blue: 3.6pm, green: 4.5 pm, red: 8.0 pm) and a Chandra adaptively smoothed false-color composite 
at right (blue: 0. 3-1.0 keV, green: 1. 0-2.0 keV, red: 2.0-10.0 keV). The green ellipses indicate the D 25 isophotes. See Sec. 2 for a brief 
description of X-ray point-source selection. 
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Figure 2. Properties of our galaxy sample compared to other SINGS galaxies. Left: Color-magnitude diagram. Right: Plot of SFR vs. 
M*. In both panels, crosses in normal font indicate SINGS galaxies that are not used in this paper. Bold crosses with diamonds indicate 
the twelve SINGS galaxies used in this paper. Galaxies are color-coded according to morphology as indicated in the legend. Values in the 
color- magnitude plot are taken from Moustakas et al. (2010), and in the SFR vs. M* plot from Noll et al. (2009). 


Table 1 

Model Parameters. 


Parameter 3, 

Notation 

Value 

Reference 

Initial Orbital Period distribution 

F(P) 

flat in log a b 

Abt (1983) 

Initial Eccentricity Distribution 

F(e) 

Thermal F(e) ~ e 

Heggie (1975) 

Binary Fraction 

fhin 

50% 


Magnetic Braking 
Metallicity 

z 

0.0001, 0.0002, 0.005, 0.001, 

Ivanova & Taam (2003) 

IMF high-end slope 


0.002, 0.005, 0.01, 0.02, 0.03 
-2.35 or -2.7 

Kroupa (2001); Kroupa & Weidner (2003) 

Initial Mass Ratio Distribution 

F(q) 

Flat, twin, or 50% flat — 50% twin 

Kobulnicky & Fryer (2007); Pinsonneault & Stanek (2006 

CE Efficiency x central concentration 

Aoce 

0.1, 0.2, 0.3, 0.5 

Podsiadlowski et al. (2003) 

Stellar wind strength 

Vwind 

0.25, 1.0, 2.0 

Belczynski et al. (2010) 

CE during HG 

CE-HG 

Yes or No 

Belczynski et al. (2007) 

SN kick for ECS/AIC C NS 


20% of normal NS kicks 

Linden et al. (2009) 

SN kick for direct collapse BH d 

«DCBH 

Yes (0.1) or No (0) 

Fragos et al. (2010) 


a The full range of values shown in column 3 is shown only for reference. Although our model grid contains results for these values, in 
this paper we only use models with the values shown in bold in column 3. As explained in the text, this corresponds to using (1) only 
models 1 — 96 (flat q distribution) and 193 — 288 (50% — 50% q distribution), and (2) only solar metallicities in all cases. 
b a is the semi-major axis of the binary orbit. 
c Electron Capture Supernova / Accretion Induced Collapse 

d The Hobbs et al. (2005) kick distribution for NS is multiplied by this parameter to introduce small kicks for BHs formed through a 
SN explosion with negligible ejected mass. 


Table 2 

SINGS galaxy subsample used in this paper. 


Galaxy 

(i) 

D 

(Mpc) 

(2) 

Hubble Type 

(3) 

T-Type 

(4) 

log SFR 
(M 0 yr- 1 ) 

(5) 

log M* 
(Mo) 
(6) 

SSFR 

(10-n yr -i) 

(7) 

-Vsrc 

(8) 

NGC1291 

10.8 

SBa 

1 

-0.66 ±0.27 

11.16 ±0.03 

0.151356 

88 

NGC2841 

14.1 

SAb 

3 

-0.31 ±0.24 

10.99 ±0.03 

0.501187 

32 

NGC3184 

11.1 

SABcd 

6 

0.11 ±0.15 

10.28 ±0.13 

6.76083 

47 

NGC3198 

13.68 

SBc 

5 

-0.03 ±0.15 

10.10 ±0.12 

7.4131 

19 

NGC3351 

9.33 

SBb 

3 

-0.04 ±0.19 

10.66 ±0.08 

1.99526 

45 

NGC3521 

10.1 

SABbc 

4 

0.43 ±0.22 

11.00 ±0.04 

2.69153 

60 

NGC3627 

9.38 

SABb 

3 

0.55 ±0.22 

10.95 ±0.05 

3.98107 

59 

NGC4631 

7.62 

SBd 

7 

0.70 ±0.20 

10.47 ±0.14 

16.9824 

29 

NGC4736 

5.2 

SAab 

2 

0.07 ±0.20 

10.80 ±0.07 

1.86209 

42 

NGC4826 

7.48 

SAab 

2 

-0.32 ±0.21 

10.77 ±0.02 

0.812831 

28 

NGC5055 

7.8 

SAbc 

4 

0.43 ±0.22 

10.96 ±0.06 

2.95121 

59 

NGC5474 

6.8 

SAcd 

6 

-0.65 ±0.14 

9.30 ±0.06 

11.2202 

16 


Note. — Columns are: (1) galaxy name; (2) distance (Moustakas et al. 2010); (3) Hubble Type 
(de Vaucouleurs et al. 1991); (4) T-type; (5) star formation rate from Noll et al. (2009); (6) loga- 
rithmic stellar mass from Noll et al. (2009); (7) specific star formation rate (SFR /M*) from (5) and 
(6); (8) number of X-ray point sources in XSINGS catalog. 
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Figure 3. Observational and theoretical cumulative XLFs for SINGS galaxies. Blue: observational (dotted, original; solid, corrected for 
incompleteness). Red: theoretical XLF that has highest pair likelihood (Table 4) for a given galaxy (dotted, original; solid, with added 
expected unrelated contribution from background AGN). Dark grey: theoretical XLF (with added expected AGN contribution) for the 
galaxy shown, but using StarTrack model 245, which has the highest global likelihood. For galaxies NGC 1291, 2841, 3351, 3627, and 4826 
the solid red and solid dark grey curves coincide. The horizontal (vertical) hashed regions indicate =blcr estimates for the oXLF and tXLF, 
respectively, due to Poisson statistics and bolometric uncertainties. 
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Figure 4. Normalized pair (top) and global (bottom) likelihood (Sec. 4.3) vs. model number for all models used in this paper. As 
explained in the text, the global curve in the bottom panel is the combination of the individual curves in the top panel. 



Figure 5. Normalized global likelihood (Sec. 4.3) vs. StarTrack parameter Aqce (donor central concentration x common envelope 
efficiency) varied in our model grid. Crosses indicate resistant mean values for normalized global likelihood values of models using a given 
Aq-ce value. Error bars indicate 3cr estimates. 
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h (erg s ') L x (erg ■*') L x (erg s ’) 

Figure 6. Cumulative theoretical XLFs for XSINGS galaxies, for subpopulations based on donor mass and age of stellar population. 
Continuous black curves show the total tXLF (black). The “old” LMXB population is shown by the red dotted lines, and the “old” HMXB 
population by the blue dotted line. The “young” LMXB population is shown by the red solid line, and the “ young” HMXB population by 
the blue solid line. The “young” and “old” subpopulations are defined in Noll et al. (2009) and are shown separately by solid and dotted 
lines, respectively (see also Table 3). Note that the Lx axes have different ranges, matching the observed XLF range. 
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DONOR TYPE (color) 

All XRBs degenerate evolved MS 


ACCRETOR TYPE (line type) 
— BH NS 






L* (erg s ') L x (erg s ') L x (erg s ') 

Figure 7 . Cumulative theoretical XLFs for XSINGS galaxies, for subpopulations based on donor and accretor type. The continuous black 
curve is the total tXLF. For colored curves, color indicates donor stellar type and line type indicates accretor type, as indicated at the top 
of the figure. 
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log SFR (M q yr -1 ) log M, (M Q ) 


Figure 8. Integrated Lx, 0 . 3 — lO.OkeV from LMXBs and HMXBs from StarTrack modeling vs. SFR and M* (Noll et al. 2009) for the 12 
SINGS galaxies in this paper. Luminosities have been calculated using both the best pair-likelihood model for each galaxy (open grey 
circles) and model 245, the best global model in this paper, which is also the top-ranked model in F13 (black filled circles). The dashed line 
shows the expected contribution from HMXBs (left panel) and LMXBs (right panel) based on the best-fit relation of L10 (Equation 10). 
The dotted lines show ±1 a uncertainties. 


Table 3 

Star formation histories from SED fitting. 


Galaxy 

(i) 

log M„ a 1 
(Mq) 
(2) 

log T oS p 
(Gyr) 
(3) 

£oSP 

(Gyr) 

(4) 

log TySP 

(Gyr) 

(5) 

*ySP 

(Gyr) 

(6) 

/burst 

(7) 

NGC1291 

11.3 

-0.60 

10.0 

-1.30 

0.20 

-2.52 

NGC2841 

11.2 

-0.60 

10.0 

-1.30 

0.20 

-2.52 

NGC3184 

10.6 

-0.60 

10.0 

0.00 

0.20 

-2.00 

NGC3198 

10.1 

1.00 

10.0 

-1.30 

0.20 

-1.52 

NGC3351 

10.9 

-0.60 

10.0 

-1.30 

0.20 

-2.00 

NGC3521 

11.2 

-0.60 

10.0 

-1.30 

0.20 

-1.52 

NGC3627 

11.1 

0.00 

10.0 

0.00 

0.20 

-2.00 

NGC4631 

10.8 

-0.60 

10.0 

-0.60 

0.20 

-1.52 

NGC4736 

11.0 

-0.60 

10.0 

0.00 

0.05 

-3.00 

NGC4826 

10.9 

-0.60 

10.0 

-1.30 

0.20 

-2.00 

NGC5055 

11.2 

-0.60 

10.0 

-1.30 

0.20 

-1.52 

NGC5474 

9.5 

0.48 

10.0 

-1.30 

0.20 

-1.00 


Note. — Best fit values from SED fitting (run B in Noll et al. (2009) 
and Noll, private communication). Columns are: (1) galaxy name; (2) 
total stellar mass plus gas mass from stellar mass loss; (3) e- folding 
timescale for old stellar population; (4) age of old stellar population; (5) 
e-folding timescale for young stellar population; (6) age of young stellar 
population; (7) mass fraction of young stellar population. 
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Table 4 

Parameters for highest pair-likelihood model-galaxy pairs 


Galaxy 

(i) 

Model 

(2) 

Aq?cb 

(3) 

IMF exponent 
(4) 

^7wind 

(5) 

CE-HG 

(6) 

q distribution 

(7) 

« DCBH 
(8) 

Rank 

(9) 

likelihood ratio 
(10) 

NGC1291 

245 

0.1 

-2.7 

1.0 

No 

50-50 

0.1 

10 

-381 

NGC2841 

245 

0.1 

-2.7 

1.0 

No 

50-50 

0.1 

1 

0 

NGC3184 

277 

0.1 

-2.7 

2.0 

Yes 

50-50 

0.1 

12 

5636 

NGC3198 

253 

0.1 

-2.7 

2.0 

No 

50-50 

0.1 

2 

-7 

NGC3351 

245 

0.1 

-2.7 

1.0 

No 

50-50 

0.1 

11 

-425 

NGC3521 

229 

0.1 

-2.7 

2.0 

Yes 

50-50 

0 

9 

-180 

NGC3627 

245 

0.1 

-2.7 

1.0 

No 

50-50 

0.1 

8 

-177 

NGC4631 

229 

0.1 

-2.7 

2.0 

Yes 

50-50 

0 

5 

-44 

NGC4736 

261 

0.1 

-2.7 

0.25 

No 

50-50 

0.1 

4 

-23 

NGC4826 

245 

0.1 

-2.7 

1.0 

No 

50-50 

0.1 

3 

-8 

NGC5055 

269 

0.1 

-2.7 

1.0 

Yes 

50-50 

0.1 

7 

-84 

NGC5474 

197 

0.1 

-2.7 

1.0 

No 

50-50 

0 

6 

-45 


Note. — Columns are: (1) galaxy ID; (2) model number; (3) CE efficiency x central concentration; (4) exponent 
of high-mass power law component of IMF: Kroupa (2001, —2.35) or Kroupa Sz Weidner (2003, —2.7); (5) stellar 
wind strength parameter; (6) Yes: All possible outcomes of a CE event with a Hertzsprung gap donor allowed; No: 
A CE with such a donor star will always result to a merger; (7) binary mass ratio distribution: 50-50 means half 
of the binaries originate in a twin binary distribution and half in a flat mass ratio distribution; (8) parameter with 
which the Hobbs et al. (2005) kick distribution is multiplied for BHs formed through a SN explosion with negligible 
ejected mass; (9) rank of best-fit galaxy-model pairs based on pair likelihood value; (10) natural logarithm of ratio 
of each pair likelihood to maximum pair likelihood in this table: ln(£ pa i r) /£pair, km, max)- 


Table 5 

Parameters and global likelihood values for 15 best models ranked by likelihood. 


Model 

(1) 

A«ce 

(2) 

IMF exponent 
(3) 

Vwind 

(4) 

CE-HG 

(5) 

q distribution 
(6) 

« DCBH 

(7) 

Rank 

(8) 

Rank F13 
(9) 

Rank T13 
(10) 

likelihood ratio 

(ii) 

245 

0.1 

-2.7 

1.0 

No 

50-50 

0.1 

1 

1 

4 

0 

253 

0.1 

-2.7 

2.0 

No 

50-50 

0.1 

2 

14 

5 

-81 

277 

0.1 

-2.7 

2.0 

Yes 

50-50 

0.1 

3 

8 

3 

-97 

229 

0.1 

-2.7 

2.0 

Yes 

50-50 

0 

4 

2 

2 

-128 

269 

0.1 

-2.7 

1.0 

Yes 

50-50 

0.1 

5 

3 

12 

-137 

205 

0.1 

-2.7 

2.0 

No 

50-50 

0 

6 

4 

1 

-198 

197 

0.1 

-2.7 

1.0 

No 

50-50 

0 

7 

13 

50 

-207 

61 

0.1 

-2.7 

2.0 

No 

Flat 

0.1 

8 

16 

7 

-258 

201 

0.1 

-2.35 

2.0 

No 

50-50 

0 

9 

10 

19 

-260 

221 

0.1 

-2.7 

1.0 

Yes 

50-50 

0 

10 

15 

60 

-262 

53 

0.1 

-2.7 

1.0 

No 

Flat 

0.1 

11 

20 

22 

-283 

273 

0.1 

-2.35 

2.0 

Yes 

50-50 

0.1 

12 

6 

6 

-285 

249 

0.1 

-2.35 

2.0 

No 

50-50 

0.1 

13 

5 

10 

-296 

85 

0.1 

-2.7 

2.0 

Yes 

Flat 

0.1 

14 

12 

9 

-305 

37 

0.1 

-2.7 

2.0 

Yes 

Flat 

0 

15 

7 

8 

-341 


Note. — - Columns are: (1) model number; (2) CE efficiency x central concentration; (3) exponent of high-mass power law 
component of IMF: Kroupa (2001, —2.35) or Kroupa & Weidner (2003, —2.7); (4) stellar wind strength parameter; (5) Yes: All 
possible outcomes of a CE event with a Hertzsprung gap donor allowed; No: A CE with such a donor star will always result to 
a merger; (6) binary mass ratio distribution: 50-50 means half of the binaries originate in a twin binary distribution and half in 
a flat mass ratio distribution; (7) parameter with which the Hobbs et al. (2005) kick distribution is multiplied for BHs formed 
through a SN explosion with negligible ejected mass; (8) rank of model based on likelihood value in this paper; (9) rank of model 
based on likelihood value in F13; (10) rank of model based on likelihood value in T13; (11) natural logarithm of ratio of each global 
likelihood to maximum global likelihood in this table: ln(£ g i 0 bai/£giobal,max)- (The full table is available on-line). 
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Table 6 

Parameters and global likelihood values for 192 models ranked by likelihood 

(on-line only). 


Model 

(1) 

Aqcb 

(2) 

IMF exponent 
(3) 

^?wind 

(4) 

CE-HG 

(5) 

q distribution 
(6) 

^DCBH 

(7) 

Rank 

(8) 

Rank F13 
(9) 

Rank T13 
(10) 

likelihood ratio 

(ii) 

245 

0.1 

-2.7 

1.0 

No 

50-50 

0.1 

1 

1 

4 

0 

253 

0.1 

-2.7 

2.0 

No 

50-50 

0.1 

2 

14 

5 

-81 

277 

0.1 

-2.7 

2.0 

Yes 

50-50 

0.1 

3 

8 

3 

-97 

229 

0.1 

-2.7 

2.0 

Yes 

50-50 

0 

4 

2 

2 

-128 

269 

0.1 

-2.7 

1.0 

Yes 

50-50 

0.1 

5 

3 

12 

-137 

205 

0.1 

-2.7 

2.0 

No 

50-50 

0 

6 

4 

1 

-198 

197 

0.1 

-2.7 

1.0 

No 

50-50 

0 

7 

13 

50 

-207 

61 

0.1 

-2.7 

2.0 

No 

Flat 

0.1 

8 

16 

7 

-258 

201 

0.1 

-2.35 

2.0 

No 

50-50 

0 

9 

10 

19 

-260 

221 

0.1 

-2.7 

1.0 

Yes 

50-50 

0 

10 

15 

60 

-262 

53 

0.1 

-2.7 

1.0 

No 

Flat 

0.1 

11 

20 

22 

-283 

273 

0.1 

-2.35 

2.0 

Yes 

50-50 

0.1 

12 

6 

6 

-285 

249 

0.1 

-2.35 

2.0 

No 

50-50 

0.1 

13 

5 

10 

-296 

85 

0.1 

-2.7 

2.0 

Yes 

Flat 

0.1 

14 

12 

9 

-305 

37 

0.1 

-2.7 

2.0 

Yes 

Flat 

0 

15 

7 

8 

-341 

225 

0.1 

-2.35 

2.0 

Yes 

50-50 

0 

16 

11 

33 

-358 

241 

0.1 

-2.35 

1.0 

No 

50-50 

0.1 

17 

31 

59 

-452 

13 

0.1 

-2.7 

2.0 

No 

Flat 

0 

18 

9 

11 

-473 

77 

0.1 

-2.7 

1.0 

Yes 

Flat 

0.1 

19 

21 

27 

-542 

261 

0.1 

-2.7 

0.25 

No 

50-50 

0.1 

20 

43 

81 

-577 

9 

0.1 

-2.35 

2.0 

No 

Flat 

0 

21 

29 

29 

-608 

57 

0.1 

-2.35 

2.0 

No 

Flat 

0.1 

22 

28 

30 

-652 

246 

0.2 

-2.7 

1.0 

No 

50-50 

0.1 

23 

25 

23 

-675 

81 

0.1 

-2.35 

2.0 

Yes 

Flat 

0.1 

24 

26 

20 

-680 

5 

0.1 

-2.7 

1.0 

No 

Flat 

0 

25 

23 

42 

-682 

29 

0.1 

-2.7 

1.0 

Yes 

Flat 

0 

26 

24 

48 

-706 

265 

0.1 

-2.35 

1.0 

Yes 

50-50 

0.1 

27 

41 

79 

-726 

254 

0.2 

-2.7 

2.0 

No 

50-50 

0.1 

28 

22 

13 

-756 

285 

0.1 

-2.7 

0.25 

Yes 

50-50 

0.1 

29 

46 

87 

-757 

206 

0.2 

-2.7 

2.0 

No 

50-50 

0 

30 

18 

14 

-757 

193 

0.1 

-2.35 

1.0 

No 

50-50 

0 

31 

56 

109 

-777 

33 

0.1 

-2.35 

2.0 

Yes 

Flat 

0 

32 

30 

32 

-782 

198 

0.2 

-2.7 

1.0 

No 

50-50 

0 

33 

39 

65 

-802 

231 

0.3 

-2.7 

2.0 

Yes 

50-50 

0 

34 

37 

31 

-831 

230 

0.2 

-2.7 

2.0 

Yes 

50-50 

0 

35 

17 

16 

-840 

278 

0.2 

-2.7 

2.0 

Yes 

50-50 

0.1 

36 

19 

15 

-872 

270 

0.2 

-2.7 

1.0 

Yes 

50-50 

0.1 

37 

34 

39 

-929 

255 

0.3 

-2.7 

2.0 

No 

50-50 

0.1 

38 

40 

21 

-933 

199 

0.3 

-2.7 

1.0 

No 

50-50 

0 

39 

55 

70 

-935 

279 

0.3 

-2.7 

2.0 

Yes 

50-50 

0.1 

40 

33 

24 

-944 

247 

0.3 

-2.7 

1.0 

No 

50-50 

0.1 

41 

45 

40 

-946 

222 

0.2 

-2.7 

1.0 

Yes 

50-50 

0 

42 

42 

71 

-990 

207 

0.3 

-2.7 

2.0 

No 

50-50 

0 

43 

32 

25 

-997 

213 

0.1 

-2.7 

0.25 

No 

50-50 

0 

44 

84 

111 

-1009 

217 

0.1 

-2.35 

1.0 

Yes 

50-50 

0 

45 

60 

115 

-1111 

49 

0.1 

-2.35 

1.0 

No 

Flat 

0.1 

46 

58 

74 

-1215 

202 

0.2 

-2.35 

2.0 

No 

50-50 

0 

47 

48 

45 

-1227 

223 

0.3 

-2.7 

1.0 

Yes 

50-50 

0 

48 

57 

85 

-1297 

274 

0.2 

-2.35 

2.0 

Yes 

50-50 

0.1 

49 

47 

41 

-1339 

262 

0.2 

-2.7 

0.25 

No 

50-50 

0.1 

50 

87 

98 

-1357 

237 

0.1 

-2.7 

0.25 

Yes 

50-50 

0 

51 

91 

113 

-1398 

208 

0.5 

-2.7 

2.0 

No 

50-50 

0 

52 

27 

18 

-1400 

73 

0.1 

-2.35 

1.0 

Yes 

Flat 

0.1 

53 

71 

88 

-1430 

226 

0.2 

-2.35 

2.0 

Yes 

50-50 

0 

54 

54 

62 

-1441 

250 

0.2 

-2.35 

2.0 

No 

50-50 

0.1 

55 

44 

36 

-1448 

271 

0.3 

-2.7 

1.0 

Yes 

50-50 

0.1 

56 

53 

56 

-1464 

93 

0.1 

-2.7 

0.25 

Yes 

Flat 

0.1 

57 

86 

107 

-1464 

69 

0.1 

-2.7 

0.25 

No 

Flat 

0.1 

58 

79 

103 

-1484 

227 

0.3 

-2.35 

2.0 

Yes 

50-50 

0 

59 

80 

73 

-1541 

232 

0.5 

-2.7 

2.0 

Yes 

50-50 

0 

60 

35 

28 

-1551 

21 

0.1 

-2.7 

0.25 

No 

Flat 

0 

61 

101 

119 

-1563 

14 

0.2 

-2.7 

2.0 

No 

Flat 

0 

62 

51 

37 

-1575 

39 

0.3 

-2.7 

2.0 

Yes 

Flat 

0 

63 

76 

57 

-1581 

62 

0.2 

-2.7 

2.0 

No 

Flat 

0.1 

64 

52 

34 

-1586 

203 

0.3 

-2.35 

2.0 

No 

50-50 

0 

65 

66 

64 

-1598 

1 

0.1 

-2.35 

1.0 

No 

Flat 

0 

66 

81 

105 

-1610 

280 

0.5 

-2.7 

2.0 

Yes 

50-50 

0.1 

67 

36 

26 

-1619 

248 

0.5 

-2.7 

1.0 

No 

50-50 

0.1 

68 

61 

55 

-1653 

256 

0.5 

-2.7 

2.0 

No 

50-50 

0.1 

69 

38 

17 

-1659 

54 

0.2 

-2.7 

1.0 

No 

Flat 

0.1 

70 

59 

58 

-1675 

38 

0.2 

-2.7 

2.0 

Yes 

Flat 

0 

71 

49 

38 

-1737 

214 

0.2 

-2.7 

0.25 

No 

50-50 

0 

72 

110 

121 

-1771 

45 

0.1 

-2.7 

0.25 

Yes 

Flat 

0 

73 

106 

122 

-1790 
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Table 6 — Continued 


Model 

(1) 

A«cb 

(2) 

IMF exponent 
(3) 

^?wind 

(4) 

CE-HG 

(5) 

q distribution 
(6) 

«DCBH 

(7) 

Rank 

(8) 

Rank F13 

(9) 

Rank T13 
(10) 

likelihood ratio 

(ii) 

25 

0.1 

-2.35 

1.0 

Yes 

Flat 

0 

74 

94 

112 

-1799 

86 

0.2 

-2.7 

2.0 

Yes 

Flat 

0.1 

75 

50 

35 

-1802 

87 

0.3 

-2.7 

2.0 

Yes 

Flat 

0.1 

76 

68 

53 

-1815 

63 

0.3 

-2.7 

2.0 

No 

Flat 

0.1 

77 

77 

46 

-1817 

257 

0.1 

-2.35 

0.25 

No 

50-50 

0.1 

78 

116 

135 

-1843 

263 

0.3 

-2.7 

0.25 

No 

50-50 

0.1 

79 

105 

114 

-1862 

6 

0.2 

-2.7 

1.0 

No 

Flat 

0 

80 

65 

72 

-1867 

251 

0.3 

-2.35 

2.0 

No 

50-50 

0.1 

81 

63 

51 

-1947 

200 

0.5 

-2.7 

1.0 

No 

50-50 

0 

82 

83 

84 

-1952 

15 

0.3 

-2.7 

2.0 

No 

Flat 

0 

83 

70 

49 

-1957 

275 

0.3 

-2.35 

2.0 

Yes 

50-50 

0.1 

84 

67 

61 

-1977 

242 

0.2 

-2.35 

1.0 

No 

50-50 

0.1 

85 

78 

83 

-2017 

7 

0.3 

-2.7 

1.0 

No 

Flat 

0 

86 

99 

89 

-2114 

78 

0.2 

-2.7 

1.0 

Yes 

Flat 

0.1 

87 

73 

67 

-2125 

55 

0.3 

-2.7 

1.0 

No 

Flat 

0.1 

88 

96 

78 

-2129 

30 

0.2 

-2.7 

1.0 

Yes 

Flat 

0 

89 

82 

82 

-2195 

224 

0.5 

-2.7 

1.0 

Yes 

50-50 

0 

90 

95 

96 

-2298 

286 

0.2 

-2.7 

0.25 

Yes 

50-50 

0.1 

91 

97 

110 

-2343 

10 

0.2 

-2.35 

2.0 

No 

Flat 

0 

92 

93 

69 

-2354 

194 

0.2 

-2.35 

1.0 

No 

50-50 

0 

93 

102 

118 

-2431 

272 

0.5 

-2.7 

1.0 

Yes 

50-50 

0.1 

94 

88 

76 

-2540 

266 

0.2 

-2.35 

1.0 

Yes 

50-50 

0.1 

95 

90 

99 

-2574 

82 

0.2 

-2.35 

2.0 

Yes 

Flat 

0.1 

96 

98 

75 

-2600 

243 

0.3 

-2.35 

1.0 

No 

50-50 

0.1 

97 

104 

97 

-2604 

204 

0.5 

-2.35 

2.0 

No 

50-50 

0 

98 

69 

63 

-2622 

215 

0.3 

-2.7 

0.25 

No 

50-50 

0 

99 

132 

130 

-2623 

35 

0.3 

-2.35 

2.0 

Yes 

Flat 

0 

100 

122 

100 

-2748 

195 

0.3 

-2.35 

1.0 

No 

50-50 

0 

101 

120 

126 

-2751 

16 

0.5 

-2.7 

2.0 

No 

Flat 

0 

102 

64 

44 

-2809 

58 

0.2 

-2.35 

2.0 

No 

Flat 

0.1 

103 

92 

68 

-2812 

31 

0.3 

-2.7 

1.0 

Yes 

Flat 

0 

104 

108 

101 

-2822 

252 

0.5 

-2.35 

2.0 

No 

50-50 

0.1 

105 

62 

47 

-2835 

34 

0.2 

-2.35 

2.0 

Yes 

Flat 

0 

106 

100 

80 

-2840 

281 

0.1 

-2.35 

0.25 

Yes 

50-50 

0.1 

107 

124 

134 

-2849 

218 

0.2 

-2.35 

1.0 

Yes 

50-50 

0 

108 

107 

124 

-2880 

11 

0.3 

-2.35 

2.0 

No 

Flat 

0 

109 

113 

93 

-2902 

276 

0.5 

-2.35 

2.0 

Yes 

50-50 

0.1 

110 

85 

66 

-2970 

40 

0.5 

-2.7 

2.0 

Yes 

Flat 

0 

111 

74 

54 

-2993 

70 

0.2 

-2.7 

0.25 

No 

Flat 

0.1 

112 

123 

131 

-3088 

238 

0.2 

-2.7 

0.25 

Yes 

50-50 

0 

113 

117 

128 

-3097 

88 

0.5 

-2.7 

2.0 

Yes 

Flat 

0.1 

114 

75 

52 

-3140 

22 

0.2 

-2.7 

0.25 

No 

Flat 

0 

115 

137 

141 

-3145 

79 

0.3 

-2.7 

1.0 

Yes 

Flat 

0.1 

116 

103 

92 

-3211 

228 

0.5 

-2.35 

2.0 

Yes 

50-50 

0 

117 

89 

77 

-3354 

64 

0.5 

-2.7 

2.0 

No 

Flat 

0.1 

118 

72 

43 

-3355 

267 

0.3 

-2.35 

1.0 

Yes 

50-50 

0.1 

119 

114 

116 

-3376 

83 

0.3 

-2.35 

2.0 

Yes 

Flat 

0.1 

120 

118 

95 

-3642 

56 

0.5 

-2.7 

1.0 

No 

Flat 

0.1 

121 

109 

94 

-3655 

287 

0.3 

-2.7 

0.25 

Yes 

50-50 

0.1 

122 

125 

129 

-3675 

59 

0.3 

-2.35 

2.0 

No 

Flat 

0.1 

123 

112 

91 

-3683 

258 

0.2 

-2.35 

0.25 

No 

50-50 

0.1 

124 

148 

153 

-3692 

244 

0.5 

-2.35 

1.0 

No 

50-50 

0.1 

125 

133 

120 

-3721 

65 

0.1 

-2.35 

0.25 

No 

Flat 

0.1 

126 

147 

155 

-3747 

196 

0.5 

-2.35 

1.0 

No 

50-50 

0 

127 

142 

138 

-3893 

239 

0.3 

-2.7 

0.25 

Yes 

50-50 

0 

128 

139 

140 

-3980 

219 

0.3 

-2.35 

1.0 

Yes 

50-50 

0 

129 

126 

133 

-3993 

209 

0.1 

-2.35 

0.25 

No 

50-50 

0 

130 

145 

149 

-3998 

264 

0.5 

-2.7 

0.25 

No 

50-50 

0.1 

131 

135 

137 

-4004 

71 

0.3 

-2.7 

0.25 

No 

Flat 

0.1 

132 

141 

145 

-4247 

89 

0.1 

-2.35 

0.25 

Yes 

Flat 

0.1 

133 

152 

156 

-4269 

8 

0.5 

-2.7 

1.0 

No 

Flat 

0 

134 

121 

106 

-4304 

50 

0.2 

-2.35 

1.0 

No 

Flat 

0.1 

135 

119 

108 

-4360 

259 

0.3 

-2.35 

0.25 

No 

50-50 

0.1 

136 

162 

166 

-4454 

94 

0.2 

-2.7 

0.25 

Yes 

Flat 

0.1 

137 

136 

139 

-4582 

216 

0.5 

-2.7 

0.25 

No 

50-50 

0 

138 

150 

144 

-4608 

210 

0.2 

-2.35 

0.25 

No 

50-50 

0 

139 

166 

160 

-4777 

2 

0.2 

-2.35 

1.0 

No 

Flat 

0 

140 

131 

127 

-4909 

23 

0.3 

-2.7 

0.25 

No 

Flat 

0 

141 

151 

152 

-4919 

32 

0.5 

-2.7 

1.0 

Yes 

Flat 

0 

142 

134 

125 

-4941 

233 

0.1 

-2.35 

0.25 

Yes 

50-50 

0 

143 

144 

148 

-4973 

74 

0.2 

-2.35 

1.0 

Yes 

Flat 

0.1 

144 

129 

123 

-5056 

12 

0.5 

-2.35 

2.0 

No 

Flat 

0 

145 

115 

90 

-5068 

46 

0.2 

-2.7 

0.25 

Yes 

Flat 

0 

146 

143 

147 

-5261 

51 

0.3 

-2.35 

1.0 

No 

Flat 

0.1 

147 

140 

132 

-5439 

26 

0.2 

-2.35 

1.0 

Yes 

Flat 

0 

148 

138 

136 

-5439 
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Table 6 — Continued 


Model 

(1) 

Aocb 

(2) 

IMF exponent 
(3) 

Vwind 

(4) 

CE-HG 

(5) 

q distribution 
(6) 

«DCBH 

(7) 

Rank 

(8) 

Rank F13 

(9) 

Rank T13 
(10) 

likelihood ratio 

(ii) 

3 

0.3 

-2.35 

1.0 

No 

Flat 

0 

149 

149 

142 

-5475 

80 

0.5 

-2.7 

1.0 

Yes 

Flat 

0.1 

150 

130 

117 

-5481 

60 

0.5 

-2.35 

2.0 

No 

Flat 

0.1 

151 

111 

86 

-5491 

84 

0.5 

-2.35 

2.0 

Yes 

Flat 

0.1 

152 

127 

102 

-5517 

41 

0.1 

-2.35 

0.25 

Yes 

Flat 

0 

153 

159 

164 

-5557 

17 

0.1 

-2.35 

0.25 

No 

Flat 

0 

154 

160 

165 

-5648 

288 

0.5 

-2.7 

0.25 

Yes 

50-50 

0.1 

155 

155 

158 

-5707 

240 

0.5 

-2.7 

0.25 

Yes 

50-50 

0 

156 

161 

163 

-5950 

268 

0.5 

-2.35 

1.0 

Yes 

50-50 

0.1 

157 

146 

143 

-6167 

220 

0.5 

-2.35 

1.0 

Yes 

50-50 

0 

158 

154 

150 

-6219 

36 

0.5 

-2.35 

2.0 

Yes 

Flat 

0 

159 

128 

104 

-6354 

211 

0.3 

-2.35 

0.25 

No 

50-50 

0 

160 

173 

170 

-6397 

282 

0.2 

-2.35 

0.25 

Yes 

50-50 

0.1 

161 

157 

157 

-6626 

75 

0.3 

-2.35 

1.0 

Yes 

Flat 

0.1 

162 

153 

146 

-6684 

18 

0.2 

-2.35 

0.25 

No 

Flat 

0 

163 

175 

176 

-7172 

47 

0.3 

-2.7 

0.25 

Yes 

Flat 

0 

164 

164 

162 

-7192 

95 

0.3 

-2.7 

0.25 

Yes 

Flat 

0.1 

165 

158 

159 

-7368 

66 

0.2 

-2.35 

0.25 

No 

Flat 

0.1 

166 

170 

174 

-7542 

27 

0.3 

-2.35 

1.0 

Yes 

Flat 

0 

167 

156 

151 

-7695 

52 

0.5 

-2.35 

1.0 

No 

Flat 

0.1 

168 

165 

154 

-7837 

4 

0.5 

-2.35 

1.0 

No 

Flat 

0 

169 

167 

161 

-7917 

234 

0.2 

-2.35 

0.25 

Yes 

50-50 

0 

170 

169 

167 

-8550 

72 

0.5 

-2.7 

0.25 

No 

Flat 

0.1 

171 

163 

168 

-8769 

24 

0.5 

-2.7 

0.25 

No 

Flat 

0 

172 

168 

169 

-9112 

283 

0.3 

-2.35 

0.25 

Yes 

50-50 

0.1 

173 

171 

172 

-9139 

67 

0.3 

-2.35 

0.25 

No 

Flat 

0.1 

174 

181 

183 

-9148 

260 

0.5 

-2.35 

0.25 

No 

50-50 

0.1 

175 

177 

178 

-9520 

212 

0.5 

-2.35 

0.25 

No 

50-50 

0 

176 

184 

181 

-10327 

19 

0.3 

-2.35 

0.25 

No 

Flat 

0 

177 

183 

184 

-10550 

90 

0.2 

-2.35 

0.25 

Yes 

Flat 

0.1 

178 

176 

177 

-11415 

48 

0.5 

-2.7 

0.25 

Yes 

Flat 

0 

179 

179 

182 

-11490 

235 

0.3 

-2.35 

0.25 

Yes 

50-50 

0 

180 

180 

175 

-11524 

96 

0.5 

-2.7 

0.25 

Yes 

Flat 

0.1 

181 

178 

180 

-11718 

28 

0.5 

-2.35 

1.0 

Yes 

Flat 

0 

182 

174 

173 

-12336 

76 

0.5 

-2.35 

1.0 

Yes 

Flat 

0.1 

183 

172 

171 

-12477 

42 

0.2 

-2.35 

0.25 

Yes 

Flat 

0 

184 

182 

179 

-12787 

284 

0.5 

-2.35 

0.25 

Yes 

50-50 

0.1 

185 

185 

187 

-15605 

236 

0.5 

-2.35 

0.25 

Yes 

50-50 

0 

186 

189 

188 

-16426 

91 

0.3 

-2.35 

0.25 

Yes 

Flat 

0.1 

187 

186 

185 

-16623 

20 

0.5 

-2.35 

0.25 

No 

Flat 

0 

188 

190 

190 

-18594 

43 

0.3 

-2.35 

0.25 

Yes 

Flat 

0 

189 

187 

186 

-19295 

68 

0.5 

-2.35 

0.25 

No 

Flat 

0.1 

190 

188 

189 

-19584 

44 

0.5 

-2.35 

0.25 

Yes 

Flat 

0 

191 

192 

192 

-29890 

92 

0.5 

-2.35 

0.25 

Yes 

Flat 

0.1 

192 

191 

191 

-30317 


Note. — Columns are: (1) model number; (2) CE efficiency x 
central concentration; (3) exponent of high-mass power law com- 
ponent of IMF: Kroupa (2001, —2.35) or Kroupa & Weidner (2003, 
—2.7); (4) stellar wind strength parameter; (5) Yes: All possible 
outcomes of a CE event with a Hertzsprung gap donor allowed; 
No: A CE with such a donor star will always result to a merger; 
(6) binary mass ratio distribution: 50-50 means half of the binaries 
originate in a twin binary distribution and half in a flat mass ratio 
distribution; (7) parameter with which the Hobbs et al. (2005) 
kick distribution is multiplied for BHs formed through a SN ex- 
plosion with negligible ejected mass; (8) rank of model based on 
likelihood value in this paper; (9) rank of model based on likeli- 
hood value in F13; (10) rank of model based on likelihood value 
inT13; (11) natural logarithm of ratio of each global likelihood to 
maximum global likelihood in this table: ln(£ g i 0 b a i/£giobal,max)- 
(The full table is available on-line). 


